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1 Abstract 

The hysteresis behaviors of anisotropic S-1 Heisenberg model have been studied 
within the effective field theory with two spin cluster. After giving the phase 
diagrams, the effect of the crystal field and anisotropy in the exchange inter¬ 
action on the hysteresis loops have been determined. One important finding 
is double hysteresis loop behavior of the system in the low temperature and 
negative crystal field region, which disappears with the decreasing anisotropy in 
the exchange interaction. This behavior was carefully investigated and physical 
explanation has also been given. 


2 Introduction 

Spin-1 (S-1) Heisenberg model has attracted some interest for many years. The 
model has exchange and single-ion anisotropies. There are many materials which 
can be handled by this type of model e.g. K2NiF/^ [T], CsNiCl^ [5], CsFeBr^ 
[3], which have weak axial anisotropy and strong planar anisotropy, respec¬ 
tively. Besides, molecular oxygen O2 adsorbed on graphite [3], NiGa2Si 0, and 
Ba^Mn20^ [6] are the more recent examples of this type of materials. Indeed, 
in real magnetic materials, single-ion anisotropy plays a major role for investi¬ 
gating the magnetic behavior of the system [7] . On the other hand, it has been 
shown from the band structure calculations that, spatially anisotropic exchange 
interaction along the x and y directions occur in several vanadium phosphate 
material systems, such as {Pb2V0{P04))2, SrZnV0{P04), BaZnVO^POi) 
and BaCdV0{P04)2 [8]. 

S-1 anisotropic or isotropic Heisenberg model has been studied by various 
methods, in order to obtain phase transition characteristics and thermodynam¬ 
ical properties. We can inspect the related literature by grouping the works as 
ground state works and finite temperature works. 

The ground-state properties of this system has been worked widely and 
quantum phase transition characteristics have been obtained. For instance, 
S-1 Heisenberg ferromagnet with an arbitrary crystal-field potential within the 
linked-cluster series expansion [^, S-1 Heisenberg antiferromagnet with uni¬ 
axial single-ion anisotropy in a field within the spin-wave approach m , two- 
dimensional quantum anisotropic S-1 Heisenberg antiferromagnet using the self 
consistent harmonic approximation m, S-1 Heisenberg antiferromagnetic spin 
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chains with exchange and single-site anisotropies in an external field within the 
density-matrix renormalization group techniques |12| . S-1 bilinear-biquadratic 
model on a honeycomb lattice with using the linear fiavor-wave (LFW) the¬ 
ory m and tensor renormalization group method m, S-l bond-alternating 
Heisenberg antiferromagnetic chain with a single-ion anisotropy in longitudinal 
and transverse magnetic fields within the infinite time evolving block decima¬ 
tion, the linearized tensor renormalization group, and the density matrix renor¬ 
malization group methods m, S-l antiferromagnetic Heisenberg model on the 
square lattice with a three-site interaction using a variety of analytical and nu¬ 
merical methods m and S-l Heisenberg model with a single-ion anisotropy on 
a triangular lattice within the the cluster mean-field approach m- 

With all these works, the nature of the ground state of the S-l Heisenberg 
model has been determined. Recently it has been shown that, the quantum 
effects can also be seen in a finite (but low enough) temperature region [T5] . 
Hence, the literature also include the works related to the quantum phase tran¬ 
sitions in the region of low temperatures, in which the quantum effects are cer¬ 
tainly still present. For example, the effect of the magnetic field on the phase di¬ 
agrams and the thermodynamical properties has been studied for a S-l bilinear- 
biquadratic Heisenberg model on the triangular lattice within the mean-field the¬ 
ory and exact diagonalizations |19j . S-l Heisenberg antiferromagnet with easy- 
axis or easy-plane single-site anisotropy on the square lattice within the series 
expansions [50], the effects of frustration between nearest, next-nearest neighbor 
and next-next-nearest neighbors of the quantum S-l anisotropic antiferromag¬ 
netic Heisenberg model on a simple cubic lattice with single ion anisotropy using 
the bond operator technique [21], are among them. Most recent studies are, the 
phase diagram of the anisotropic S-l Heisenberg chain with single ion anisotropy 
(D) using a ground-state fidelity approach [22] . the low-temperature properties 
of one-dimensional S-l Heisenberg model with geometric fluctuations within the 
strong-disorder renormalization-group and quantum Monte Carlo and density- 
matrix renormalization-group numerical calculations [23] and two-dimensional 
S-l antiferromagnet with next nearest neighbor exchange interactions and easy 
axis single ion anisotropy, on the square lattice, are studied at low temperature 
using a Modified Spin Wave Theory [51] . 

On the other hand, when we look at the phase diagrams and thermodynam¬ 
ical properties of the system in the finite temperature region, we can see several 
works with various methods, as in the ground properties of the system. 

One of the first studies in this case is, S-l uniaxial ferromagnetic model with 
both exchange anisotropy and single-ion anisotropy (H), as well as transverse 
coupling are studied in the mean-field approximation [25] . The phase diagrams 
of the system have been obtained in this work. One-dimensional S-l ferromag¬ 
netic isotropic Heisenberg model is studied by the double-time Green’s function 
method [26] and two-dimensional classical isotropic Heisenberg antiferromagnet 
with a single-ion anisotropy is studied in the presence of a uniform magnetic field 
along the easy axis within the MC simulation and Green’s function technique 

m- 

The three components of the magnetizations for the anisotropic Heisenberg 
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model with single ion anisotropy have been calculated by the use of many-body 
Green’s function method, when external field is applied in both x and z di¬ 
rections [35]. The effect of an easy-plane crystalline anisotropy and easy-axis 
exchange anisotropy on the phase diagram of the three-dimensional classical 
ferromagnetic anisotropic Heisenberg model has been studied by using Monte 
Carlo simulations |29j . In addition, the thermodynamic properties of S-1 ferro¬ 
magnetic chains with an easy-axis single-ion anisotropy have been investigated 
by both a Green-function approach, based on a decoupling of three-spin opera¬ 
tor products, and by exact diagonalizations of chains with up to iV = 12 sites 
using periodic boundary conditions [50] . 

The general works related to spin S Heisenberg systems have also been 
present. The thermodynamical properties of one and two-dimensional ferro- 
magnets with arbitrary spin-S in a magnetic field have been investigated by a 
second-order Green-function theory for the isotropical Heisenberg model |31j . 
Also the phase diagrams of the anisotropic ferromagnetic spin-S Heisenberg 
model on a square lattice have been obtained by double-time Green’s function 
method within the Callen decoupling approximation |32j . 

Most recent studies related to the thermodynamic properties and phase tran¬ 
sition characteristics are, S-1 Heisenberg antiferromagnet with easy-plane single¬ 
ion anisotropy on three-dimensional bipartite lattices with sixth-order series 
expansions |33j . three-dimensional anisotropic Heisenberg XXZ model with a 
crystal field by using the variational approach based on the Bogoliubov inequal¬ 
ity [34] and quantum S-1 anisotropic ferromagnetic Heisenberg model within 
the same method [35j . The field-induced laws of thermodynamic properties are 
obtained by Green’s function method for the two-dimensional S-1 ferromagnetic 
Heisenberg model with the exchange and single-ion anisotropies [36j . Some ex¬ 
act solutions were aslo reported such as S-1 Ising-Heisenberg diamond chain in 
a magnetic held by a rigorous treatment based on the transfer-matrix method 

ISZl- 

As seen in this short literature, the problem of thermodynamical properties 
of the S-1 anisotropic Heisenberg model is up to date. Thus, the aim of this 
work is to determine the hysteresis behaviors of the S-1 anisotropic Heisenberg 
model on a bulk system. For this aim, the paper is organized as follows: In Sec. 
|3|we briefly present the model and formulation. The results and discussions are 
presented in Sec. H] and hnally Sec. [5] contains our conclusions. 

3 Model and Formulation 

We start with a standard S-1 anisotropic Heisenberg Hamiltonian as, 

n = -jJ2 -^E*^" ( 1 ) 

<i,j> i i 

where Sf,Sf,Sf denote the x,y,z component of the Pauli S-1 operator at a 
site i respectively, J stands for the exchange interactions between the nearest 
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neighbor spins, A is the anisotropy in the exchange interaction, D and H are 
the crystal field and longitudinal external magnetic field at sites of the lattice, 
respectively. The first summation is carried over the nearest neighbors of the 
lattice, while the others are over all the lattice sites. We note that, Hamiltonian 
given in Eq. (HD covers Ising model, anisotropic Heisenberg model and isotropic 
Heisenberg model depending on the value of A. When A = 0, Eq. ([T]) represents 
the S-1 Ising model (Blume-Capel model) while A = 1 gives the isotropic S-1 
Heisenberg model. The intermediate values of the anisotropy in the exchange 
interaction (0 < A < 1) corresponds to the anisotropic S-1 Heisenberg model, 
by means of the XXZ model. 

In an EFT-2 approximation [5S], we start by constructing the 2-site cluster 
and write 2-site cluster Hamiltonian with the axial approximation |39j as 

2 2 

= -J [A + SfS^^) + 5^51] - dY, iSD" - E (2) 

i=l i=l 

Here hi denotes all the interactions between the spin at a site i and the other 
spins which are outside of the cluster, where i = 1, 2. Let the site i has number 
of nearest neighbor sites which are located at the outside of the selected 
cluster, then hi can be written within the axial approximation [39) as 


h, = jYsIu. 


k^l 


(3) 


where is the k. the nearest neighbor of the spin Sf. The thermal average 
of the quantity S'f (i = 1,2) via exact generalized Callen-Suzuki identity [JD] is 
given by 


/ Trzgf exp(-/3H(")) \ 

\ rr 2 exp (—y 


(4) 


In Eq. (U), Tr 2 stands for the partial trace over the lattice sites 1 and 2, which 
belong to the constructed finite cluster and /3 = l/(fc_BT), where ks is the 
Boltzmann constant and T is the temperature. Calculation of Eq. requires 
the matrix representation of the related operators in chosen basis set, which can 
be denoted by where i = 1,2,..., 9. Each of the element of this basis set 

can be represented by |siS 2 ), where Si = ±1,0 is just one spin eigenvalues of 
the S-1 operator Sf {i = 1,2). 

Let us denote the matrix elements of the operator defined in Eq. 

by, 

Rij — (^4^ i 

These elements can be found in Sec. El 

In order to get the matrix representation of the operator exp (—the 
matrix whose elements defined in Eq. (ED has to be diagonalized. But the diag¬ 
onal form of this matrix cannot be obtained analytically, thus some numerical 
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procedures have to be applied. After numerical diagonalization, diagonal ele¬ 
ments of this matrix (where a = 1, 2,..., 9) can be obtained. Needless to 
say, these are just the eigenvalues of the matrix whose elements are dehned in 
Eq. (O. Let us denote the transformed basis set as |(/)') (i = 1, 2,..., 9) which 
makes the matrix representation of the operator diagonal, 


7 ^( 2 ) 


( 6 ) 


Let us denote the diagonal matrix elements of the operator S'f in this new 
basis set as 


= (7) 

With these definitions, Eq. O can be written as 

9 

a exp {—j3ra) , 

{St) = l—, -y (8) 

^ exp(-/3r„) 

The order parameter of the system can be defined by 

^ = \{{Sl) + {S^2))- ( 9 ) 

By writing Eq. ([8|) in Eq. we can get the magnetization of the system 
in a closed form as 

m = {f {hi,h2)), (10) 

where the function is defined as 
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{Tl^a + 72,a) exp (-/3ra) 

fihuh2) = l^ -^-. ( 11 ) 

exp i-PVa) 

We note that the function in Eq. (ED also the functions of the Hamiltonian 
variables (J, A, H, D) as well as the temperature. Needless to remind that, the 
function in Eq. m cannot given analytically. 

Evaluating Eq. (nni with using differential operator technique and decou¬ 
pling approximation |41j is possible. This approximation is most widely used 
for these type of systems within the effective field theory. Within this technique 
Eq. m can be written as 

m = (exp(/iiVi + h2V2)) fixi,X2)\xi=o,x2=o ( 12 ) 
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where Vi = V 2 = are the differential operators, and the effect of the 
exponential differential operator on an arbitrary function /(xi, X 2 ) is defined by 

exp (aiVi + a 2 V 2 )/(xi, X 2 ) = /(xi + oi, X 2 + 03 ), (13) 

where oi, 02 are the arbitrary constants. 

With writing Eq. ([3|) in Eq. (TT^ and using S-1 Van der Waerden identity 
|42j . we can arrive the expression of the order parameter as 


m = nn [1 + sinh (JVj) + (cosh (JV,) - 1)] \ /(xi, X 2 )Ui=o, 2 : 2=0 

/ 

(14) 

The quadrupolar moment of the system which is defined by 

1 ^ (((5'^)^) + <^(S'|)^)) = (exp(/iiVi + h2V2))5(xi,X2)|xi=o.x2=o (15) 

can be calculated in the same way of the magnetization of the system as 


q = 


^ Zj, 

nn 




[1 + S'i.fc sinh (JVj) + Sf i^ (cosh (JVi) - 1)] 


where the function defined by 


, ^ 2 ) laii—0,ai2 —0 
(16) 
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{Ul,a + U 2 ,a) exp {-Pro) 

g{hi,h 2 ) = ]~ -g-. (17) 

exp(-/3ra) 


Here [7^,^ is given by 


Ui a — 



(18) 


We note that by writing Eqs. (ITTl) and (fTO . we have assumed that the spins Si 
and S '2 have no common nearest neighbors. 

If we use the decoupling approximation [41] in Eqs. o and (fTHI) and assume 
the homogeneity of the lattice (i.e. all lattice sites are equivalent) we get the 
expressions 


m = 0f 02 "/(a;i,X 2 )Ui=o,x 2 =o 
q ~ 01^ 02^ff (^1: ^2) |xi=0,a:2=0 

where the operator is defined by 

0i = [1 + msinh {JVi) + q (cosh {JVi) — 1)]. 


(19) 


( 20 ) 
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Eq. m is a coupled nonlinear equation system of m, q. Before the numerical 
solution, last step contains writing the operators in Eq. dm as exponentials in 
order to use Eq. m- Hence, Eq. dm can be written as 




r. 

n—Op—0 

with Binomial expansion. Here the coefficients are defined by 


C^l = 


Zl 

n 

r+p 


n—p r p 

EEE 

r=0 


n—p 

r 


/ \ r+p 

f-j (-1)" P ’'+''exp[(r+p-2t-2r;) JVi]. 


( 21 ) 


( 22 ) 


The critical temperature of the system can be obtained by linearizing Eq. 
(fTOl) in m. 
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4 Results and Discussion 


In this section, after the phase diagrams of the system is reviewed, hysteresis 
characteristics of the system will be obtained, for a square lattice (zi = Z 2 = 3) 
as an example of the 2D lattices. We use the scaled quantities given as 


keT 
J ’ 



D 

7 


(23) 


within the calculations. 

4.1 Phase diagrams 

The phase diagrams of the square lattice within the S-1 anisotropic Heisenberg 
model can be seen in Fig. [1] for selected values of anisotropy in the exchange 
interaction (A). We can see from Fig. [T]that, all curves related to the different 
anisotropy in the exchange interactions have tricritical points, at which plotted 
phase diagrams end. At this point, the second order magnetic phase transition 
meets the first order one. We note that, when the anisotropy in the exchange 
interaction decreases (i.e. A increases), the tricritical point shifts towards the 
low temperature and crystal field region of the (tc, d) plane. We can say that 
all curves have the same characteristics, regardless of the value of A : when 
the crystal field decreases from the positive region, the critical temperature 
of the system decreases. This decreasing behavior of the critical temperature 
of the system becomes more rapid in the negative crystal field region. At a 
specific value of the crystal field, second order critical curves terminated (at 
specific tricritical values). The isotropical model has lower critical temperature 
than the anisotropical model, for all values of d. Since rising anisotropy in the 
exchange interaction (decreasing A) forces the spins align along the z axis, then 
it is not surprising that, the required energy supplied by the temperature to 
destroy the order of the spins z direction, has to be higher. 

We can compare the results in the critical temperatures with other methods 
for some known limits of the model such as d —> 00 (namely Ising limit), d = 0. 
The critical temperature values obtained in the present work can be seen in 
Table [U in comparison with the other methods. We can see from Table [1] that, 
the critical temperatures lie between the variational approach |35| and mean 
field approximation (MFA) [25]. But we can say that the limiting results of 
Ising model can be considered as good values in its own class, namely in EFT. 
The critical value of S-I/2 limit (i.e. d —>• 00) is tc = 3.025 is slightly better than 
the tc = 3.090 which is most widely used formulation namely DA m within 
the EFT. Besides, when we compare the value of the critical temperature in the 
case d = 0, A = 0, we can see that this value is slightly lower than the value 
obtained within the DA |4l], which is tc = 2.187. Thus it is expected that, the 
critical values are slightly lower than the DA [41] for the anisotropic Heisenberg 
model. Nevertheless we have to say that, the result of the isotropical Heisenberg 
limit, contradicts to the Mermin-Wagner Theorem [45] . according to which it is 
not possible to observe any long range order in isotropic Heisenberg model in 
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Figure 1: The variation of the critical temperature with crystal field for selected 
values of A, for anisotropic S-1 Heisenberg model on a square lattice. 


Table 1: Critical temperatures in some limiting cases for the square lattice in 
comparison with corresponding exact, variational and MFA values. 



Exact 

Variational [55] 

Present 

MFA [ID 

d —>■ oo 

2.269[43] 

2.885 

3.025 

4.0 

d = 0,A = 0 

1.693[4| 

2.065 

2.149 

2.667 

d = 0,A = l 

0.0[45] 

1.492 

1.722 

2.667 


one and two dimensions. But, this deficiency in limiting case can not prevent 
the investigation of the effect of the anisotropy in the exchange interaction on 
the hysteresis loop behaviors, which is the main topic of this work. 

4.2 Hysteresis Behaviors 

In order to investigate the behavior of the hysteresis loops with the variation 
of the Hamiltonian parameters, we can use quantities which are related to the 
shapes of the hysteresis loops, namely hysteresis loop area (HLA), coersive field 
(CF) and the remanent magnetization (RM). Since the evolution of the hystere¬ 
sis loops with the variation of the temperature is trivial, we inspect the effect 
of varying crsytal field (d) and the anisotropy in the exchange interaction A on 
the hysteresis loops. 

Increasing temperature enhances the thermal fluctuations and this causes to 
destruct the magnetic order of the system. Hence the HLA, CF and RM values 
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tend to vanish and ferromagnetic hysteresis loops evolve into paramagnetic loops 
by increasing temperature. 

Let us remind the well known effect of the crystal held (d) on the order 
parameter of the system. With given value of A, in case of d > 0, system will 
order along the easy axis z up to critical temperature value tc- But for the 
values that provide d < 0, the z axis becomes hard direction and when the 
value of d lowers, the spins tend to align in xy plane, which can be regarded as 
easy plane. Thus we can say that, the disordered phase of the system can be 
related to the random alignment of the spins in 2 direction (d > 0 and t > tc), 
as well as random alignment of spins in xy plane (d < 0). 

In order to see the effect of the crystal held on the hysteresis behavior, 
we depict the variation of the HLA with crystal held, for selected values of 
anisotropy in the exchange interaction in Fig. [2l for selected temperatures as 
(a) t = 0.1 and (b) t = 0.5. The same plots for the CF and RM can be seen in 
Figs. Oandm respectively. 




Figure 2: The variation of the hysteresis loop area with crystal held, for selected 
values of A, at temperature (a) t = 0.1 and (b) t = 0.5, for anisotropic S-1 
Heisenberg model for square lattice. 

In Fig. ([2]), the curves labeled by A,B,C correspond to the S-1 Ising model, 
XXZ model and S-1 isotropic Heisenberg model, respectively. For positive values 
of the d, varying A has almost no effect on the HLA. At hrst sight, a qualitative 
distinction between the curves A and C stand out in relief. The constant HLA 
takes place in case of Ising model at large negative values of the crystal held, 
while this is not the case in isotropic Heisenberg model. By inspecting the 
corresponding CF and RM curves (see the curves labeled by A in Figs. ([3]) (a) 
and (ID) (a) at large negative values of crystal held) we can say that this nonzero 
HLA do not originate from the symmetric hysteresis loops around the origin. 
Because of the large negative values of the crystal held, the system cannot stay 
in the ferromagnetic phase, then the hysteresis loops of this region correspond 
to the paramagnetic phase. Since corresponding CF and RM values are zero at 
large negative values of the crystal held, with nonzero HLA there has to be loops 
which are not symmetric about the origin in the m — h plane. The hysteresis 
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loops for selected Hamiltonian parameter values from this large negative valued 
crystal field region can be seen in Fig. [5] 




Figure 3: The variation of the coersive field with crystal field, for selected values 
of A, at temperature (a) t = 0.1 and (b) t = 0.5, for anisotropic S-1 Heisenberg 
model for square lattice. 




Figure 4: The variation of the remanent magnetization with crystal field, for 
selected values of A, at temperature (a) t = 0.1 and (b) t = 0.5, for anisotropic 
S-1 Heisenberg model for square lattice. 

In Fig. [SJ we depict the hysteresis loops calculated at the temperature 
t = 0.1 for crystal field values (a) d = 4.0, (b) d = —1.0, (c) d = —2.2 and (d) 
d = —4.0 for a square lattice. Each plot contains two representative examples 
of hysteresis loops corresponding to A = 0.0 (Ising model, solid black curve) 
and A = 0.5 (anisotropic Heisenberg model, dotted red curve). First of all, we 
see from Fig. [S]that, the loops related to the anisotropic Heisenberg model lie 
inside of the loops for the Ising model. This is consistent with the results shown 
in Figs. [2]|4l which conclude that, decreasing anisotropy (i.e. rising A) causes to 
a decline in the HLA, RM and CF. The remarkable result seen in Fig. [5] (d) is 
double hysteresis behavior. When the crystal field decreases, first the hysteresis 
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loops get narrower (compare loops in Figs. [5] (a) and (b) ), then after a specific 
value of d, loops start to split up (see[5](c)), after this double hysteresis loop 
behavior survives. This result can explain the constant and nonzero behavior 
of the hysteresis loops with decreasing d, at constant temperature. 

Since the large negative value of d forces the spins to align in xy plane, the 
double hysteresis loop behavior, which appears in large negative values of the 
crystal field may come from the forced alignment of the spins from this plane 
to the z direction. We can see this situation by defining a quantity (1 — g), 
which is the measure of the number of spins aligned in xy plane. In Fig. [5] we 
depict the ’’hysteresis” loops for the 1 — g, with same parameters used in Fig. 
[5l We can see from Figs. [5] (a) and (b) that 1 — g of the system always zero, 
regardless of the value of the magnetic field. All spins are aligned along parallel 
to z direction. But in contrast to this situation, as seen in Fig. [6](d), magnetic 
field can induce the transition to alignment of the spins from the xy plane to 
the z axis, when the value of the magnetic field large enough. This situation is 
consistent with the results given in Fig. [5] 

We can summarize these results as: the decreasing crystal field gives rise to 
double hysteresis behavior. This double hysteresis loop behavior comes from 
the transition between the states which are consist of aligned spins in xy plane 
and spins in z direction. This mechanism is different from the mechanism which 
explains the hysteresis behavior in positive valued crystal field. 

The double hysteresis loops have been observed experimentally in different 
systems, for example in Mn- doped (pb, La) {Zr,Ti)Oz ceramics [4^, and in 
FezOi/Mn^O^ super lattices [47]. In addition, triple hysteresis behaviors have 
been observed experimentally, such as single chain magnets with antiferromag¬ 
netic interchain coupling [48] and in molecular-based magnetic materials [49], 
CoFeB/Cu, CoNip/Cu, FeGa/py, and FeGa/CoFeB multilayered nanowires 

m- 

In theoretical explanations, this double (or triple) hysteresis loop behaviors 
mostly attributed to the exchange interaction ratios in the nano material. For 
instance, the effect of the transverse field on the hysteresis behavior of the S-1 
Ising nanotube has been investigated within the EFT based on a probability 
distribution method. Double and triple hysteresis loops were obtained m- 
Another work was devoted to the hysteresis behavior of the S-1 Ising ferrielectric 
cubic nanowire with negative core shell coupling within the MC simulation. It 
is observed that, when the absolute value of J shell-core / J core increases, then the 
hysteresis curve changes from one central loop to triple loops [52] ■ Also, it was 
shown within the EFT formulation based on a probability distribution method 
that, trimodal random field distribution on the S-1 Ising nanotube could give 
rise to double hysteresis loops when the system passes from the ferromagnetic 
phase to a paramagnetic one [S3]. In addition to such studies, we show in 
the present work that, double hysteresis loop behavior can also occur in bulk 
systems as explained above. Lastly, one important study deals with the S- 
llsing model with transverse crystal field (in x direction) within the EFT with 
probability distribution technique |54| . In this work the authors cannot find 
double hysteresis loop behavior due to the transverse alignment of the crystal 
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field. 

Lastly, the double hysteresis loop region can be obtained and it is plotted 
in Fig. [71 In Fig. [71 closed loop seperate the single and double hysteresis loop 
regions in {t, A) plane. It can be seen from Fig. jTjthat, double hysteresis loops 
may occur in low temperature and high anisotropy in the exchange interaction 
region. After a value of A = 0.591/t = 0.630, system cannot show double 
hysteresis loop behavior at any t/A. In other words, decreasing anisotropy in 
the exchange interaction (rising A) or rising temperature destroys the double 
hysteresis loop behavior. 
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Figure 5: Hysteresis loops of the XXZ model with A = 0.5 (shown in dotted 
red) and Ising model (shown in solid black) on square lattice with temperature 
t = 0.1 for crystal field values (a) d = 4.0, (b) d = —1.0, (c) d = —2.2 and (d) 
d= -4.0. 


5 Conclusion 

The effects of the crystal field and anisotropy in the exchange interaction of 
the anisotropic S-1 Heisenberg model on the hysteresis behavior have been ob¬ 
tained within the effective field theory with two spin cluster. The model has 
three Hamiltonian parameters as crystal field (d), anisotropy in the exchange 
interaction (A) and the temperature (t). The model covers the Ising model 
(A = 0.0 or d —oo) and isotropic Heisenberg model (A = 1.0) as limiting 
cases. 

After the phase diagram of the system for several values of A given, the 
effects of d and A on the hysteresis loops have been determined. The value of 
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Figure 6: Hysteresis loops of the XXZ model with A = 0.5 (shown in dotted 
red) and Ising model (shown in solid black) on square lattice with temperature 
t = 0.1 for crystal field values (a) d = 4.0, (b) d = —1.0, (c) d = —2.2 and (d) 
d= -4.0. 

A makes no important difference in the hysteresis loops in low t and positive d 
region. However in the negative d region one fundamental distinction appears 
for hysteresis loops that have different values of A, namely double hysteresis loop 
behavior. Double hysteresis loop behavior appears for low values of t and A, 
when d takes large negative values. This behavior is related to the alignment 
of spins in xy plane and tendency of this alignment to the 2 direction when 
external longitudinal magnetic field is applied. It is shown that, this behavior 
cannot appear in the isotropic Heisenberg model, indeed for the values that 
provide A > 0.591. Double hysteresis loop region depicted in (t, A) plane has 
also been investigated in detail and physical explanation given briefly. 

We hope that the results obtained in this work may be beneficial form both 
theoretical and experimental point of view. 
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Figure 7: Double and single hysteresis loop regions of the S-1 anisotropic Heisen¬ 
berg model on a sqnare lattice in t — A plane. The cnrve obtained for a large 
negative value of the crystal field. 

Appendix: Matrix elements of the Hamiltonian 

Matrix elements of the 2-spin anisotropic Heisenberg Hamiltonian given by Eq. 
m was defined by Eq. ([S]). The nonzero matrix elements are as follows: 


i?ii = -J-{2D + 2H + hi+h2) 
R22 = -(D + iJ + hi) 

i?24 ~ —2TA 

i?33 = J — {2D hi — ^ 2 ) 

i?35 = —2JA 
i?42 = —2TA 

i?44 = — {D H -\- h2) 

R53 = —2JA 
i?57 = —2JA 
Ree = -{D-H-h2) 

Res = —2JA 
i?75 = —2TA 
i?77 = J — {2D — hi -\- / 12 ) 

Rse = —2JA 

Rss = -{D-H-hi) 

Rgg = —J — {2D — hi — / 12 ) 
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